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ABSTRACT 

Aims. We seek to estimate the average level of MHD turbulence near coronal mass ejection (CME) fronts as they propagate from the 
Sun to the Earth. 

Methods. We examined the cosmic ray data from the GRAPES-3 tracking muon telescope at Ooty, together with the data from other 
sources for three closely observed Forbush decrease events. Each of these event is associated with frontside halo coronal mass ejections 
(CMEs) and near-Earth magnetic clouds. The associated Forbush decreases are therefore expected to have significant contributions 
from the cosmic-ray depressions inside the CMEs/ejecta. In each case, we estimate the magnitude of the Forbush decrease using 
a simple model for the dilfusion of high-energy protons through the largely closed field lines enclosing the CME as it expands 
and propagates from the Sun to the Earth. The diffusion of high-energy protons is inhibited by the smooth, large-scale magnetic field 
enclosing the CME and aided by the turbulent fluctuations near the CME front. We use estimates of the cross-field diffusion coefficient 
Dj^ derived from the published results of extensive Monte Carlo simulations of cosmic rays propagating through turbulent magnetic 
fields. We then compare our estimates with the magnitudes of the observed Forbush decreases. 

Results. Our method helps constrain the ratio of energy density in the turbulent magnetic fields to that in the mean magnetic fields 
near the CME fronts. This ratio is found to be ~ 2% for the 2001 April 1 1 Forbush decrease event, ~ 6% for the 2003 November 20 
Forbush decrease event and ~ 249% for the much more energetic event of 2003 October 29. 

Key words, cosmic rays; Sun: coronal mass ejection (CME); solar-terrestrial relations 



1. Introduction 

. Forbush decreases are short-term depressions in the cosmic ray 
' flux reaching the Earth, and they are caused by the effects of the 
interplanetary counterparts of coronal mass ejections (CMEs) 
from the Sun (and the shocks they drive) and also the co- 
rotating interaction regions originating from the Sun. They have 
, been studied closely since their initial discovery in the 1930s 
' (see, e.g., comprehensive observational reviews by Cane 2000; 
Venkatesan & Badruddin 1990). With the recent upsurge of in- 
terest in space weather effects due to solar transients, the com- 
plementary information provided by the cosmic ray signatures 
of these effects has assumed increased significance. 

In this work, we have examined the data from the GRAPES- 
3 tracking muon telescope at Ooty for three well-observed 
events. Earlier, we had analysed the Forbush decrease data from 
this experiment during the period 2001-2004 and searched for 
events that can be associated with a near-Earth magnetic cloud 
and a corresponding halo CME from the Sun. We start from 
lists of magnetic cloud events given on the WIND webpage, 
http://lepmfi.gsfc.nasa.gov/mfi/magxloud_publ.html, Lynch et 
al. (2003) and Huttunen et al. (2005). We shoitlist Forbush de- 
crease events that occurred in a reasonable time window around 
the time of entry of the magnetic clouds. Of these shortlisted 
events, we have selected the three events shown in Fig. |8] Two 
of these events (2001 April 11 and 2003 October 29) have es- 
pecially clean decrease profiles. Even though the decrease pro- 
file for the third event on 2003 November 20 is not nearly as 



clear-cut, we still selected it, because of its association with a 
very closely observed magnetic cloud, and also because it illus- 
trates the diversity of the events observed with the GRAPES-3 
experiment. We estimate the contribution of the depressed cos- 
mic ray density inside the CME to the total Forbush decrease 
seen in these datasets, and use it to derive conclusions regarding 
the level of turbulence near the CME front. 

The events discussed in this work were very closely ob- 
served, and they have their origins in full halo CMEs that orig- 
inated close to the centre of the solar disc. This means that 
the corresponding interplanetary CMEs (ICMEs) would have 
been intercepted at the Earth as ejecta/magnetic clouds, and 
the observed cosmic ray depressions coiTesponding to these 
events would have had contributions from the shock ahead of 
the ICME, as well as from the ejecta/magnetic cloud itself 
(Cane et al. 1994; Cane et al. 1995; Cane & Richardson 1997). 
Theoretical treatments of Forbush decreases model the effect as 
arising due to a general propagating region of enhanced tur- 
bulence/scattering and decreased diffusion (e.g., Nishida 1983; 
le Roux & Potgieter 1991) and do not distinguish between 
the shock and the ejecta, or implicitly assume that the de- 
crease is only due to the shock (e.g., Chih & Lee 1986). It is 
fairly well known that magnetic clouds are well-correlated with 
Forbush decreases (e.g., Zhang & Burlaga 1988; Badruddin et 
al. 1986; Badruddin et al. 1991; Venkatesan & Badruddin 1990; 
Sanderson et al. 1990; see, however, Lockwood et al. 1991 for 
the opposite viewpoint). 
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The relationship between the decrease due to the shock and 
the one due to the magnetic cloud has been investigated (Cane et 
al. 1995, 1997; Wibberenz et al. 1997, 1998), and there is some 
evidence that they can contribute in roughly equal proportions, 
to the overall magnitude of the decrease, although the individual 
time profiles of the decrease due to these two effects can be quite 
different (Wibberenz et al 1997). It is worth mentioning that the 
contributions of the shock and the magnetic cloud to the overall 
decrease can be discerned relatively easily from the spacecraft 
observations (e.g.. Cane et al. 1995, 1997; Cane 2000), since the 
temporal boundaries of the magnetic cloud can be readily ascer- 
tained by using the data from the same source. These data typi- 
cally exhibit two steps in the decrease phase, one corresponding 
to the shock passage and the second to the entry of the magnetic 
cloud. The two-step signature, however, need not be a neces- 
sary condition for the Forbush decreases arising from the com- 
bination of a magnetic cloud and the associated shock. Even for 
events where there is clear evidence that both the shock and the 
following magnetic cloud have intercepted the Earth (as is the 
case with the three events presented in this paper), the precise 
decrease profile will depend upon the standoff distance between 
the magnetic cloud and the shock and on the energies of the pro- 
tons detected. If the protons are energetic enough that their mean 
free path is comparable to the magnetic cloud-shock standoff 
distance, the distinct identities of the steps corresponding to the 
shock and the cloud may not be clear. We attribute 50% of the 
total decrease to the magnetic cloud (a = 0.5, see § 3), follow- 
ing the general logic outlined by Wibberenz et al. (1997). We 
discuss the sensitivity of our final results to variations in a in 
§ 6. In our quantitative work, we estimate the magnitude of the 
Forbush decrease using filtered data (see § 2.1 for a description 
of the filtering method) averaged over 1 hr. 

The high muon counting rate measured by the GRAPES-3 
experiment results in extremely small statistical errors, allow- 
ing small changes in the intensity of the cosmic ray flux to be 
measured with high precision. Thus a small drop (~0.2%) in the 
cosmic ray flux, during a Forbush decrease event can be reli- 
ably detected. This is possible even in the presence of the di- 
urnal anisotropy of much larger magnitude (~1.0%), through a 
suitable filtering technique described subsequently. 

Before we proceed further, it is worth discussing the ter- 
minology we use. The term CME is generally used to denote 
the blob of plasma ejected from the solar corona as viewed by 
coronographs near the Sun (typically from 1 to 30 Rq). The CME 
can be thought of as a largely closed, magnetic flux rope-like 
magnetic structure (see, however, Bothmer et al. 1996, 1997 for 
evidence of CMEs possessing closed, as well as some open field 
lines) whose cross-section expands as it travels through the he- 
liosphere. There are practically no high-energy cosmic rays in- 
side it when it starts out from near the Sun, and the ambient cos- 
mic rays difiiise into it via cross-field diff'usion across the closed 
magnetic field lines as it travels through the heliosphere. The 
term "ICME" is a broad one used to denote the interplanetary 
counterpart of a CME (see, e.g., Forsyth et al. 2006; Wimmer- 
Schweingruber et al. 2006 and references therein), while the 
term "magnetic cloud" is reserved for ICMEs detected near the 
Earth, possessing certain well-defined criteria such as plasma 
temperature depressions and smooth magnetic field rotations 
(see, e.g., Burlaga et al. 1981, Bothmer & Schwenn 1998). Other 
CME/ICME related ejecta detected near the Earth are merely 
called "ejecta". 

The rest of this paper is organised as follows. We briefly de- 
scribe the GRAPES-3 experiment in § 2. We describe the model 
we used to calculate the magnitude of the expected Forbush de- 
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Fig. 1. A schematic layout the GRAPES-3 air shower array 
shown here with 217 detectors (circles). Each of the 16 squares 
represent a 35 m^ area tracking muon detector with E/j > 1 GeV 
used in the present work. 



creases in § 3. The three Forbush decrease events we analysed 
are described in § 4, 5, and 6. We describe our results in § 7 and 
summarise them in § 8. 



2. GRAPES-3 Experimental system 

The GRAPES-3 experiment is located at Ooty (1 1.4°N latitude, 
76.7°E longitude, and 2200 m altitude), in southern India. The 
GRAPES-3 air shower experiment was designed to have a com- 
pact configuration of a conventional type array, with a separa- 
tion of only 8 m between the adjacent detectors, which are de- 
ployed in a symmetric hexagonal geometry. A schematic layout 
of the GRAPES-3 array is shown in Fig. [1] The observations 
were started in 2000 with 217 detectors, located within the inner 
8 rings, see Fig.[T](Gupta et al. 2005). 

A very large-area tracking muon telescope operating as a 
part of the GRAPES-3 experiment (Gupta et al. 2005; Hayashi 
et al. 2005), is a unique instrument, to search for the high-energy 
protons emitted during the active phase of a solar flare or a 
coronal mass ejection (CME). The muon telescope is capable 
of providing a high-statistics, directional study of muons. The 
GRAPES-3 muon telescope covers an area of 560 m^, consist- 
ing of a total 16 modules, each 35 m^ in area. These modules are 
located close to each other as shown in Fig. [T] A cluster of four 
35 m^ area neighbouring modules, located inside a common hall, 
constitutes one supermodule with a total area of 140 m^. The en- 
ergy threshold of the telescope is 1 GeV for the muons arriving 
along the vertical direction. The cutoff' rigidity due to the mag- 
netic field of the Earth at Ooty is 17 GV in the vertical direction 
and varies from 12 to 42 GV across the field of view of the muon 
telescope as shown later in Fig.|6] 

The basic detector element of the muon telescope is a rugged 
proportional counter (PRC) made from a 600 cm long steel pipe 
with 2.3 mm wall thickness and a square cross-sectional area of 
10 X 10 cm^. A muon telescope module with a sensitive area of 
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Fig. 2. A schematic display of the 4-layer tracking muon tele- 
scope module with 58 PRCs per layer. The four layers of the 
PRCs labelled Layer-0, Layer- 1, etc. are embedded in concrete 
blocks. Inchned lines represent a set of parallel muon tracks. 



35 m^ consists of a total of 232 PRCs arranged in four layers of 
58 PRCs each, with alternate layers placed in orthogonal direc- 
tions. Two successive layers of the PRCs are separated by a 15 
cm thick concrete layer, consisting of 60x60x15 cm^ blocks as 
shown in Fig.|2] The four-layer PRC configuration of the muon 
modules allows a 3-D reconstruction of the muon track direction 
to an accuracy of ~ 6°. It is worth noting that the accuracy grad- 
ually increases with increasing zenith angle due to the greater 
separation of the triggered PRCs. 

To achieve an energy threshold of 1 GeV for vertical muons, 
an absorber of total thickness ~ 550 g cm"^ in the form of con- 
crete blocks is employed. This was done by placing a total of 15 
layers of concrete blocks above Layer- 1 , as shown schematically 
in Fig. |2] A unique feature of the GRAPES-3 muon module is 
the robust structure of the PRCs, which supports the huge load 
of 2.4 m thick concrete in a self-supporting manner. The con- 
crete blocks are arranged in the shape of an inverted pyramid 
to shield the PRCs, with coverage up to 45° around the vertical 
direction for the incident muons. The threshold energy changes 
to sec 9 GeV for the muons incident at a zenith angle of 0. The 
cross section of a muon telescope module is shown schemati- 
cally in Fig. 121 A cluster of four such modules, separated by a 
horizontal distance of 130 cm at the base constitutes one super- 
module. The GRAPES-3 muon telescope contains a total of four 
supermodules (Hayashi et al. 2005). 

Because of the sensitivity of the PRCs to low-energy j rays, 
from the radioactivity present in the concrete absorber, individ- 
ual PRCs display sizable counting rates of ~ 200 Hz. An output 
is generated, if any one of the 58 PRCs produces a signal. A 
logical 'OR' of outputs from all 58 PRCs in a layer is gener- 
ated, after suitable amplification and shaping to form the layer 
OR output. A coincidence of the four OR outputs from the four 
layers in a module is used to generate the 4-layer trigger. Despite 
high counting rates of individual PRCs, the 4-layer coincidence 
trigger is relatively free of the contribution from the background 
radioactivity and it is caused only by the passage of a muon. The 
observed 4-layer muon counting rate of ~ 3200 Hz per module 
yields a total counting rate ~ 3 x 10^ min"' for all 16 modules. 
This high rate permits small changes of < 0.1% in the muon 
flux, to be detected on a timescale of ~ 5 min, after appropriate 
correction to the variation in the atmospheric pressure with time. 

Most of the detected muons are generated by > 20 GeV 
galactic cosmic rays, and form a stable and dominant back- 
ground to the variation in their flux produced by the CME/solar 
flare. The muon data is grouped online every 10 sec, into solid- 
angle bins of ~ 0.05 sr, consistent with the angular resolution of 
the muon telescope as described below. Since > 1 GeV muons 
are secondaries produced by the primary protons of energy > 20 
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Fig. 3. Distribution of the relative variation in %, for the 4-layer 
muon counting rate for the 35 m^ area modules for time intervals 
of 6 minutes. Also shown is a Gaussian fit to the data. The stan- 
dard deviation of the fit is 0.103% as against an expected value 
of 0.095% based on statistics. 



GeV in the atmosphere, therefore these observations can be used 
to probe the eff'ect of the Sun on cosmic rays. 

The data recorded during 2003 October, from 15 out of 16 
working modules were used in the statistical analysis described 
below. The muon rate for each 35 m^ area module was recorded 
for 770 time intervals, each 6 minutes in duration spread over a 
total of about 3.5 days. During this period of 3.5 days, the data 
could not be recorded for 70 out of 840 intervals, due to the 
failure of the recording system in the first supermodule. Next, 
the mean muon rate of the module for the entire duration was 
computed. Then the percent deviation from the mean module 
rate was calculated, separately for each module. Finally, the rms 
spread in the percent deviation called 'relative variation' was cal- 
culated for each module for all 6-minute intervals. In Fig.|3] the 
distribution of the relative variation in the 4-layer muon rate in 
percent, as recorded by the muon modules is shown. Also shown 
in Fig. |3] is a Gaussian fit to this distribution, which yields a 
standard deviation of 0. 103%. This value may be compared with 
an expectation of 0.095%, based on a 6-minute muon statistic 
of 1.1 X 10* per module. This tiny statistical error in muon rate 
allows a high-precision study of various solar phenomena to be 
carried out (Nonaka et al. 2006). 



2.1. Data analysis 

In all four supermodules, the direction of muons is recorded 
into 225 solid-angle bins. This was done by using a dedicated 
direction-sensitive trigger with an independent data acquisition 
system for each of the four supermodules. As shown in Fig. H] 
the muon angle is determined for each PRC in the lower layer 
and binned into 15 angular bins based on the specific location of 
the PRC triggered in the upper layer from among the 15 PRCs, 
one directly above (central PRC) and 7 each on either side of the 
central PRC. This angular binning is carried out in each of the 
two orthogonal projections (XZ and YZ; Z is vertical direction), 
thereby generating a 2-dimensional 15 x 15 - 225 solid-angle 
map of muon directions. The contents of the 225 solid-angle 
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Fig. 4. A schematic view of muon arrival angle selection based 
on the PRC triggered in the lower and 15 PRCs in the upper 
layer. The triggered PRCs are shown as filled squares. 
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Fig. 5. A schematic of the nine muon arrival direction bins; 3x3 
vertical bin V, and four 3x5 central bins N, E, W, S, and four 
5x5 outer bins NE, SE, SW, NW. 



bins are recorded, once every 10 sec, thus providing a contin- 
uous monitoring of the directional flux of muons in the sky. 

The variation in the muon rate may be studied in any of the 
225 solid angle bins. However, it is expected that the influence 
of a solar flare and/or CME would be spread over several bins. 
This directional spread could arise from the influence of the ter- 
restrial, solar, and interplanetary magnetic fields, etc. Therefore, 
the detected muons have been regrouped into 3x3=9 coarse 
solid-angle bins, as shown schematically in Fig.|5] This regroup- 
ing of the data was done by combining either a set of 3 x5 or 5 x5 
fine solid-angle bins. The exception being the vertical direction 
where 3x3 bins have been combined. This choice of angular 
segmentation was dictated by the fact that the muon flux is com- 
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Fig. 6. The 9 coarse solid-angle bins are shown along with the 
contours of constant geomagnetic cutoff rigidity in the field of 
view (FOV). Cutoff' rigidity varies from 12 to 42 GV in the FOV 
of GRAPES-3. 



paratively larger for the near central directions (N, E, W, S) than 
for the outer directions (NE, SE, NW, S W). This choice results in 
a relatively similar solid-angle coverage for the nine coarse bins. 
Thus the solid-angle of acceptance includes only 13x13 = 169 
out of original 225 bins, restricting the maximum zenith angle 
to 50°. This still exceeds the shielding coverage of 45°, for the 
PRCs at the outer edge, but such events constitute <1% of the 
data. This regrouping also results in muon statistics for various 
bins that are not too dissimilar. 

At the energies (> 20 GeV) of interest here, the propagation 
of charged particles near the Earth (< 20Re, where Re is the 
radius of the Earth) is strongly influenced by the geomagnetic 
field. The access by a charged particle to a given geographical 
location depends on the momentum per unit charge of the parti- 
cle called rigidity. The threshold value of the rigidity is termed 
"geomagnetic cutoff rigidity", which depends on the geograph- 
ical location on the Earth and the direction of the arriving par- 
ticle. The geomagnetic cutoff rigidity can be calculated using 
a detailed model of the geomagnetic field (Cooke et al. 1991). 
The geomagnetic cutoff rigidity for the field of view (FOV) of 
the GRAPES-3 muon telescope varies significantly for the nine 
coarse, solid-angle bins. We first calculate the cutoff rigidity for 
the centre of each of the 169 fine solid-angle bins, which con- 
stitute the 9 elements of the FOV, using the IGRF2000 geomag- 
netic field model. Subsequently, a weighted mean of the cutoff 
rigidities of the fine bins constituting a coarse bin is calculated 
for each of the 9 coarse bins. These weights are the muon count- 
ing rates for a given fine bin. Here, it needs to be emphasised that 
the knowledge of the geomagnetic field at any particular moment 
in time is imperfect, and it is virtually impossible to determine 
the cutoff to a high degree of accuracy. However, the calculated 
values of the geomagnetic cutoff represent a very useful approxi- 
mation to the true values at the time of the observations. In Fig.|6] 
the contours of constant geomagnetic cutoff rigidity in the FOV 
are superimposed over a schematic of the 9 solid-angle bins of 
muon arrival directions. The geomagnetic cutoff rigidity varies 
from 12 GV in the west to 42 GV in the east, wifliin the FOV of 
the GRAPES-3 tracking muon telescope. 
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Frequency (day 
Fig. 7. The power spectrum of 28 day data covering the period 
from 26 March to 22 April 2001 . The solid lines show the results 
for all 9 directions used in the study while the dashed line shows 
the function that is used for filtering out the high-frequency com- 
ponents. The Fourier transform is multiplied by this function 
before taking the inverse Fourier transform to get filtered time 
series. 



In this study we used the GRAPES-3 data summed over a 
time interval of one hour for each of the nine bins, which are 
identified as NW, N, NE, W, V, E, SW, S, and SE. The calcu- 
lated cutoff rigidities for these nine bins are 15.5, 18.7, 24.0, 
14.3, 17.2, 22.4, 14.4, 17.6, and 22.4 GV, respectively as seen 
from Fig. |6] The summing over an interval of 1 hour improves 
the signal-to-noise ratio, but the diurnal variations in the muon 
flux are still present. To filter out these oscillations we applied 
a low-pass filter, removing all frequencies higher than 1 day 
In Fig. |7] we show the Fourier transform of the data covering 
a period of 28 days that includes the event of 2001 April 11. 
The peaks corresponding to diurnal variation and its first har- 
monic are visible. The Fourier transform is multiplied by the 
function shown by a dashed line to remove the high-frequency 
components. This filter is found to be effective in removing high- 
frequency oscillations, including the diurnal variations and their 
harmonics. The Forbush decrease events are also clearer in the 
filtered data. 

Although the smoothing may tend to change the amplitude 
of the decrease and possibly shift the onset time for the Forbush 
decrease by a few hours in some cases. However, it is often dif- 
ficult to determine if these differences are artifacts of smoothing 
or whether the unfiltered data showed different amplitude be- 
cause a diurnal oscillation happened to have the right phase, so 
as to enhance or reduce the amplitude of the Forbush decrease. 
Some fluctuations in muon flux could be due to Forbush decrease 
and associated events, but it is unlikely that these will be peri- 
odic in nature, hence are not likely to be affected by the filter. 
Since the differences in amplitude caused by filter are not sub- 
stantial, it would not affect our results significantly. Thus in this 
work, we used the filtered data to study the characteristics of 
Forbush decrease events, although the unfiltered data are also 



shown in corresponding figures for direct comparison. Further, 
we repeated the calculations using unfiltered data and find that 
the final results are not significantly different. There can also be 
anisotropics intrinsic to the CME/magnetic cloud itself, arising 
from aBxV N drift, where B is the interplanetary magnetic field 
and V A' denotes the cosmic-ray density gradient inside the CME 
(Bieber & Evenson 1998; Munakata et al. 2003, 2005; Kuwabara 
et al. 2004). Such anisotropics can potentially be "mixed" with 
the diurnal anisotropy, and it is possible that there will still be 
some residual anisotropy after the filter is applied. This is the 
case with one of the weaker events that we have studied here; 
namely 2003 November 20. 

To determine the properties of Forbush decrease, we divided 
the time interval in three parts, the first part before the onset; the 
second part includes the decrease, while the third part includes 
the recovery phase. We fit straight lines to the first two parts and 
an exponential profile for the third part. The point of intersection 
of the two straight lines defines the onset time for the Forbush 
decrease, while the intersection of the last two fits determines 
the minimum of the profile. The difference in the flux between 
these two points determines the magnitude M of the decrease 
(§ 3), which is measured in term of percentage change in steady 
flux. The fits for all events used in this work are shown in the 
respective sections. 

In the subsequent analysis we discuss in detail three 
Forbush decrease events, namely those on 2001 April 11, 2003 
October 29, and 2003 November 20, which were associated 
with CME/magnetic clouds. The GRAPES-3 data is recorded 
along nine directions as described earlier However, in Fig. |8] 
the data for only the vertical direction is shown for these 
three events along with measurements from the Tibet neutron 
monitor at Yangbajing (Muraki et al. 2007) as obtained from 
http://crO.izmiran.rssi.ru/tibt/main.htm The Tibet neutron mon- 
itor has a wide solid-angle of acceptance with the dominant con- 
tribution coming from the vertical direction; for that reason, we 
selected the GRAPES-3 data along the vertical direction for a 
comparison as shown in Fig. [8] The Tibet neutron monitor was 
chosen, as Yangbajing (90.5°E) is located at a similar longitude 
as Ooty (76.7°E), and the geomagnetic cutoff rigidities are also 
not too dissimilar (Ooty=17 GV, and Yangbajing=14 GV). The 
neutron monitor data was scaled down by a factor of 3 to fit the 
same scale. Here, it needs to be emphasised that both sets of data 
were subjected to the same low pass filter to remove frequencies 
higher than 1 day The detailed features in the Tibet neutron 
monitor data follow the GRAPES-3 muon data rather closely for 
all three Forbush decrease events as seen from Fig. [8] 

3. Proton diffusion into CIVIE/magnetic cloud 

We are concerned here with the component of the Forbush de- 
crease produced by the difference between the proton density 
inside and outside the magnetic cloud that intercepts the Earth. 
To do this, we obtain an estimate of the proton density in- 
side the CME/magnetic cloud produced by the cumulative ef- 
fect of protons diffusing into the closed magnetic structure of 
the CME/magnetic cloud as it propagates towards the Earth. We 
use cross-field diffusion coefficients derived from recent numeri- 
cal treatments of cosmic ray propagation through turbulent mag- 
netic fields. 

Before doing so, it is worth taking a brief look at the typical 
proton gyro-radius in relation to the overall size of the ejecta to 
see how far such a picture involving cross-field diffusion would 
be valid. This is especially pertinent here, since the protons de- 
tected by the GRAPES-3 experiment are typically very ener- 
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Fig. 8. Data from Tibet neutron monitor (dashed lines) is compared with GRAPES-3 data (solid line) for vertical direction for the 
3 events studied in this work. Both the time series are filtered by applying a low pass filter to remove high frequency components. 
In order to fit in the same scale the neutron monitor data is scaled down by a factor of 3. 



getic. The GRAPES-3 experiment can detect muons produced 
by protons of rigidities above 12-42 GV, where the rigidity Rg 
is defined by 



Pc 

Rg (volts) = — = 300 B (gauss) (cm) , 
Ze 



(1) 



where P is the proton momentum, c the speed of light, Z the 
charge state (= 1 for a proton), e the charge of an electron, 
B the magnetic field in Gauss, and r/. the proton gyro-radius 
in cm. The magnetic field of a typical near-Earth magnetic 
cloud is ~ 10 G, while its radius is ~ 0.2 AU, or ~ 3 x 
lO'^ cm (jhttp://lepmfi.gsfc.nasa.gov/mfi/magj;loud_Avg.html||. 
The gyro-radius of a 30 GV proton in the field of a typical near- 
Earth magnetic cloud is ~ 10'^ cm (from Eq[T]), or 0.3 times 
the magnetic cloud radius. Estimates of the magnetic field of 
a CME near the Sun are notoriously hard to come by, but the 
rough lower bound of ~ 0. 1 G given by Bastian et al. (2001) can 
be taken as a working number. The gyro-radius of a 30 GV pro- 
ton in such a magnetic field is ~ 3 x 10*^ cm, which is around 
3 orders of magnitude smaller than the typical size of a CME at 
such heights, which is a few solar radii (a few times 10" cm). 
Thus the gyro-radius of the typical high-energy proton whose 
signature is detected by the GRAPES-3 array ranges from 10"^ 
to 0.3 times the macroscopic size of the structure into which it 
penetrates, and the picture of cross-field diffusion is generally 
valid. 

The flux of protons entering the CME/magnetic cloud at a 
given time is 



F(cm-2s-i) 

or 



(2) 



where F is in units of number per cm^ per sec, is the cross- 
field diffusion coefficient, and Ng the ambient density of high- 
energy protons. This diffusion takes place throughout the cross- 
section of the expanding CME, so the total number of high- 
energy protons that will have diffused into the CME after a time 
Tis 



A{t)F{t)dt= D^A{t)—^dt, 
Jo or 



where A{t) is the cross-sectional area of the CME at time t. The 
integration extends from the time the CME is first observed in 
the LASCO field of view {t - 0) through the time {t - T) when 
it arrives at the Earth as a magnetic cloud. A reasonable estimate 
for dNa/dr would be 



dr ~ R{t) ' 



(4) 



where R{t) is the CME cross-sectional radius at time t. 
Assuming the CME to be an expanding "flux-rope" whose 
length increases with time as its cross-sectional area expands 
(see http://lepmfi.gsfc.nasa.gov/mfi/Mag_Cloud_Model.html for 
a cartoon), we can use 



A(f) = 2nR{t)L(f), 



(5) 



where L{t) is the length of the flux-rope cylinder at time f, and 
can be related to the "height" H of the CME above the solar limb 



L(f) = 2nH(t). 



(6) 



The height above the solar limb H is the most commonly 
measured quantity for limb CMEs, but it is not clear how well 
it can be estimated for the halo CMEs we are considering. We 
consider this aspect further in § 4 below. Using Eqs. (|5]), Q, and 
^ in Eq. ([3]), we have 



'a r 

Jo 



Hit) dt . 



The density inside the near-Earth magnetic cloud would be 



nR{TfL(T) ' 



(7) 



(8) 



and it can be related to the magnitude M of the Forbush decrease 
via 



(3) a;M = 



Na-Ni AN , An D^H(f)dt 



Na Na 



= 1 - 



R^-(T) L(T) 



(9) 
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where L{T) and RiT) denote the length and cross-sectional ra- 
dius at time T, when the ejecta has reached the Earth. The quan- 
tity a denotes the fraction of the total decrease that can be at- 
tributed to the CME/magnetic cloud. As mentioned earlier, this 
fraction is typically around 50% (Wibberenz et al. 1997), so we 
adopt a = 0.5 in our calculations. We discuss the effect of varia- 
tions in a on our final result in § 7. 



3. 1 . Cross-field diffusion coefficient 

We now turn our attention to the cross-field diffusion coefficient 
D^, which governs the diffusion of the ambient high-energy pro- 
tons into the CME across the magnetic fields that enclose it. 
Protons can travel perpendicular to the direction of the average 
magnetic field as a result of normal transport of particles due to 
scattering or drift, as well as the braiding/mixing of the mag- 
netic field itself. The extent of cross-field diffusion depends on 
several parameters, such as the proton rigidity (which indicates 
how tightly the proton is bound to the magnetic field) and the 
level of turbulence, which can contribute to field line transport, 
as well as proton scattering. 

The issue of the cross-field diffusion coefficient is a subject 
of intensive ongoing research, and a good way to proceed for our 
purposes is to use the convenient parametrisation of the results 
of extensive Monte Carlo simulations of cosmic rays propagat- 
ing through tangled magnetic fields given by Candia & Roulet 
(2004). Their results can be regarded as a superset of similar 
work (e.g., Giacalone & Jokipii 1999; Casse et al. 2002); it con- 
tains all the preceding results and also extends the parameter 
regime of the simulations. They give the following fit for the 
"parallel" diffusion coefficient Dy (which is due to scattering of 
the particles back and forth along the mean field, as the field is 
subject to random turbulent fluctuations): 



(10) 



where c is the speed of light and the quantities N\\, y py 
are constants specific to different kinds of turbulence whose 
values are given in Table 1 of Candia & Roulet (2004). The 
quantity L^ax is the maximum length scale of the turbulence. 
Matthaeus et al. (2005) quote a value of Lmax = 10^ km, based 
on multi-spacecraft measurements of solar wind turbulence near 
the Earth. This is similar to values quoted by other authors (e.g., 
Manoharan et al. 1994) who use very different methods. We take 
the quantity L^ax to be a fixed number from near the Sun to 
the Earth. Candia & Roulet (2004) call p the "rigidity", but it 
is somewhat different from our usual definition. It is related to 
the usual rigidity by 

rL Rg 




(11) 



where Bq is the mean magnetic field. The quantity cr^ is the tur- 
bulence level, and is defined as 



(12) 



where is the turbulent magnetic field and the angular braces 
denote an ensemble average. The cross-field diffusion coefficient 
{D^} is related to the parallel one {D\\) by 



p < 0.2 
p > 0.2 . 



(13) 



The quantities A^^ and are constants specific to different 
kinds of turbulent spectra, and are given in Table 1 of Candia 
& Roulet (2004). For concreteness, we assume a Kolmogorov 
turbulence spectrum in our calculations. Equations ( fTOb and ( [13]) 
jointly define the cross-field diffusion coefficient. The magnetic 
field Bo of the CME changes (weakens) as it propagates from 
the Sun to the Earth, which means that p varies with time. We do 
not know how Bq varies with heliocentric distance, but we can 
assume that the flux is frozen into the CME as it propagates from 
the Sun to the Earth. It may be argued that some of the magnetic 
flux will be dissipated in driving the CME. However, on average, 
even the flux detected by near-Earth magnetic clouds is enough 
to account for some 74% of what is required to drive the CME 
from the Sun to the Earth (Subramanian & Vourlidas 2007). The 
frozen-flux assumption is therefore fairly accurate. This means 
that the product of the magnetic field and the cross-sectional area 
of the CME remains approximately constant, i.e.. 



Bo(t) - Bmc 



R(T) 



R(t) 



(14) 



where Bmc is the magnetic field inside the magnetic cloud. The 
quantity R{t) is related to the measured expansion velocity Vexp 
and the starting radius Rq of the halo CME by 



R{t)^V,,pt + Ro 



(15) 



We note that the measured quantities Rq and R{t) are the lateral 
extents of the halo CME measured in the plane of the sky, and 
are thus representative of the lateral expansion of the CME cross- 
section. 

3.2. Constraining o-^ 

For a given Forbush decrease event, we measure the quantity 
M, and Eqs. ([Hi, 0, and ^ jointly define the 
cross-field diffusion coefficient D^, which can be used in Eq. (|9]l. 
Unfortunately, there is not much simplification possible in the 
expression for D^; the quantity p typically straddles regimes 
such that both the terms under the square root in Eq. ( fTOl i and 
both the branches of Eq. (fT3T l need to be retained, and when D± is 
finally substituted into Eq. (|9]l the integral needs to be evaluated 
numerically. As we see below, we have concrete observational 
values for all the quantities used in Eq. (|9]l except the turbulence 
level cr- (EqfTSTi. and it is this quantity that we seek to constrain. 
Owing to the different cutoff rigidities in different bins and re- 
maining anisotropics, the magnitude M of the Forbush decrease 
(for a given event) differs from bin to bin. Consequently, we ob- 
tain slightly different estimates for cr^ for a given event, which 
can be regarded as a scatter (Tables 1, 2, and 3). For a given bin, 
we adopt a single number for the proton rigidity, which is equal 
to the cutoff rigidity for that bin. For instance, we adopt Rg = 17 
GV for the vertical bin, and the values of Rg for the other bins 
correspond to the respective cutoff rigidities listed in § 2.1. We 
note that this is a simplifying assumption, for a muon detected 
in a given bin can in fact be produced by protons possessing 
a range of energies. We next describe the characteristics of the 
three events we have selected for our study. 



4. 2001 April 11 event 

4.1. Magnetic cioud 

This event had its genesis in a halo CME that started at 05:30 
UT on 2001 April 10 ( ,http://cdaw.gsfc.nasa.gov/CMEiist/J . The 
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Table 1. Derived parameters for Forbush decrease for the 2001 April 1 1 event 



Quantity 


NW 


N 


NE 


W 


V 


E 


SW 


S 


SE 


FD magnitude 


2.51 


2.51 


1.87 


2.92 


2.71 


2.04 


2.58 


2.37 


1.79 


FD onset 


15:36 


17:31 


19:26 


10:05 


11:45 


16:19 


06:58 


09:22 


12:45 


FD end' 


15:07 


13:55 


12:43 


12:57 


12:43 


11:45 


09:36 


11:02 


11:17 


(Kolmogorov) 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 



' The FD end times refer to 2001 April 12. 



starting radius was Rq - 2.84 /?0, and the lateral (plane-of- 
sky) expansion speed was Vexp = 2411 kms"'. Based on an 
examination of several halo CMEs, Schwenn et al. (2005) con- 
clude that the radial speed V is typically 0.88 times the expan- 
sion speed Vexp of a halo CME. This gives V - 0.88 Vexp - 
2121 km s"'. The corresponding magnetic cloud was detected 
by the ACE spacecraft at 04:00 UT on 2001 April 12, while 
the shock reached the Earth at 15:18 UT on 2001 April 11 
(Lynch et al. 2003). Munakata et al. (2003, 2005) have mod- 
elled the magnetic cloud associated with this event as a cylinder 
with a Gaussian cross-sectional profile for the cosmic -ray den- 
sity depletion inside it. They analysed the cosmic -ray directional 
anisotropy from muon detector measurements from several sta- 
tions. Assuming that this anisotropy arises from a diamagnetic 
drift due to the interplanetary magnetic field and the density gra- 
dient inside the magnetic cloud, they have derived best-fit pa- 
rameters for the 3D cylindrical model of the magnetic cloud as 
it passes across the Earth. The length of the flux rope at time t is 

L{t)^2nH{t)^2n{Vt + Ri)). (16) 

The length L at time T is simply 

L(r) = 27rlAU. (17) 

It may be noted that we are using the value 7?o for the initial 
value of the observed height H above the solar limb. This is not 
quite right for full halo CMEs, since Rq is really the observed 
lateral extent of the CME. However, this is the only concrete 
observational quantity available, and the error is not expected 
to contribute significantly to our final result for cr^ (see § 7). 
The time elapsed between the first observation of the halo CME 
and the detection of the magnetic cloud by the ACE spacecraft 
is 46.5 hours. The average speed of the magnetic cloud by the 
time it reaches the ACE spacecraft is 640.6 kms"' (Lynch et 
al. 2003), and the spacecraft is located around 1.5 x 10^ km from 
the Earth. We therefore estimate that the magnetic cloud would 
have taken ~ 0.6 hours to traverse the distance between ACE and 
the Earth. The total time T — 46.5 H-0. 6 — Al A hours. The radius 
of the magnetic cloud is R{T) - 0.106 AU, and magnetic field 
at the centre of the cloud is estimated to be (Lynch et al. 2003) 
Bmc = 17.8 nT. 

4.2. Forbush decrease 

Table 1 gives the start and end times and magnitudes of the 
Forbush decrease in different bins for this event, as well as the 
computed values of the turbulence level cr^ (EqfT^li. The shock 
associated with this event impacted the Earth at 15:18 UT, 2001 
April 1 1 (Huttunen et al. 2005) and the magnetic cloud first in- 
tercepted the ACE spacecraft at 04:00 UT, 2001 April 12 (Lynch 
et al. 2003). Evidently, the magnetic cloud was observed on the 
Earth after the start of the Forbush decrease and well before its 
end in all the directional bins. It therefore stands to reason that 
it will have contributed substantially to the overall Forbush de- 
crease. As discussed in § 3, we adopt a value of c = 0.5, which 




10 15 20 10 15 20 10 15 20 



Fig. 9. The muon flux in the nine directions is shown for the 
Forbush decrease on 2001 April 11. The fluxes are shown as 
percentage deviation from mean values. The dashed red lines 
show the unfiltered data, while solid black lines show the data 
after applying a low-pass filter The dotted blue lines show the 
fits to filtered data. The vertical lines in each panel mark the 
arrival times of shock and magnetic cloud at the Earth. 



means that the magnetic cloud contributes to 50% of the overall 
Forbush decrease. 

5. 2003 October 29 event 

5.1. Magnetic cloud 

This was an exceptional event associated with an X17 flare 
close to the central meridian, together with a very fast halo 
CME. This full halo CME was first observed at a starting 
radius of Rq = 5.84 /?o at 11:30 UT on 2003 October 28 
(http://cdaw.gsfc.nasa.gov/CMEJist/). As with the other events, 
the quantity 7?o really represents the first observed lateral extent 
of the halo CME. The quantity L{T) can be taken to be the same 
as defined in Eq. ( fTTT i. The lateral expansion speed of the CME 
as discerned from coronograph data is Vexp - 2459 kms ' , and 
the radial expansion speed is related to Vexp by V = 0.88 Vexp as 
before. The length of the flux rope at a given time t is given as 
before by Eq. ( fTSl l. 

Kuwabara et al. (2004) have derived best-fit parameters for 
the 3D cylindrical model of the magnetic cloud as it passes 
across the Earth using the methods described in Munakata et 
al. (2003, 2005). They estimate the radius of the magnetic cloud 
to be R{T) - 0.14 AU, the magnetic field at the centre of the 
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Quantity 


NW 


N 


NE 


W 


V 


E 


SW 


S 


SE 


FD magnitude 


6.81 


6.78 


5.21 


7.76 


7.76 


5.85 


6.92 


6.82 


5.26 


FD onset 


00:58 


23:46* 


22:34* 


00:58 


00:00 


23:17* 


00:53 


00:14 


23:31* 


FDend 


18:00 


17:16 


16:19 


19:12 


18:14 


16:19 


18:43 


17:16 


15:07 


cr- (Kolmogorov) 


2.48 


2.49 


2.55 


2.45 


2.45 


2.52 


2.48 


2.49 


2.55 



* These times refer to 2003 October 28. 
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Fig. 10. The muon flux in each of the nine directions for the 
Forbush decrease on 2003 October 29. The hne-styles and mark- 
ings are the same as Fig.|9] 



where the magnitude of Forbush decrease derived here from the 
fit to the fihered data are systematically lower than those ob- 
served in the raw data, which were reported earlier (Nonaka et 
al. 2006). The onset and recovery times in the 2003 October 29 
FD event were relatively faster and thus the reduction in mag- 
nitude is largely caused by the removal of a higher frequency 
component due to the use of a low-pass filter to remove the diur- 
nal variations. The shock ahead of the ejecta impacted the Earth 
at 05:58 UT, 2003 October 29 (Dryer et al. 2004), and the mag- 
netic cloud was detected by the ACE spacecraft at 12:00 UT, 
2003 October 29. We noted earlier that it takes around 0.3 hours 
for the cloud to travel the distance between ACE and the Earth at 
a speed of 1400 km s"'; in addition to this, if one accounts for the 
travel time between the leading edge of the cloud and the centre 
using the same speed (and a cloud radius of 0.14 AU), it results 
in an additional time of 4 hours and 21 minutes. When added 
to 12:00 UT, the final result agrees very well with the intercept 
time of 16:27 UT quoted by Kuwabara et al. (2004) for the cen- 
tre of the cloud. As with the 2001 April 11 event, the magnetic 
cloud intercepts the Earth after the start of the Forbush decrease 
and well before its end in all the directional bins. The density 
depletion inside the magnetic cloud will therefore account for 
a substantial part of the overall Forbush decrease. This justifies 
our use of a - 0.5 (§ 3) to account for the fraction of the overall 
Forbush decrease related to the CME/magnetic cloud. 



cloud to be Bmc = 44 nT , and the time of closest approach of 
the cloud axis to be 16:27 UT on 2003 October 29. The ACE ob- 
servations (Huttunen et al. 2005) reveal that the magnetic cloud 
made its entry at 12:00 UT, 2003 October 29. The time elapsed 
between the first observation of the halo CME (1 1:30 UT, 2003 
October 28) and the time it was intercepted at the ACE space- 
craft (12:00 UT, 2003 October 29) is 24.5 hours . 

We could not find data for the magnetic cloud speed as mea- 
sured by the ACE spacecraft. We therefore used the estimate of 
~ 1400 km s"' from Kuwabara et al. (2004) for the magnetic 
cloud speed. At this speed, the magnetic cloud would take only 
~ 0.3 hours to traverse the distance of 1.5 x 10^ km between 
the ACE spacecraft and the Earth. The total time T is therefore 
T ^ 24.8 hours. It may be noted that the CME Hftoff speeds 
for this event and one previously discussed (2001 April 11) are 
similar However, the event of 2001 April 11 decelerates signif- 
icantly en route to the Earth, as a result of which the total travel 
time T is almost twice as long as for this event. 

5.2. Forbush decrease 

The Forbush decrease on 2003 October 29 was among the largest 
observed by the GRAPES-3 muon telescope, and it is described 
in detail in Nonaka et al. (2006). The time profiles of the de- 
crease as observed by the different modules is shown in Fig. 10 
of that paper. The derived parameters from the Forbush decrease 
profiles in the different directional bins are given in Table 2, 



6. 2003 November 20 event 

6. 1 . Magnetic cloud data 

This magnetic cloud event was associated with an erupting fil- 
ament and a halo CME, and was observed very well (Wang et 
al. 2006; Huttunen et al. 2005) by the WIND and ACE space- 
craft. Wang et al. (2006) have modelled this cloud as a cylin- 
drical, force-free flux rope to determine its radius and other pa- 
rameters. As we shall see, the Forbush decrease associated with 
this event is less than those for the other two events studied in 
this paper However, it illustrates the diversity of events observed 
with the GRAPES-3 experiment, and is associated with the best- 
studied magnetic cloud of the three events considered in this pa- 
per 

The radius of the magnetic cloud as measured by the ACE 
spacecraft is (Wang et al. 2006) RiT) = T (7 hours) x V,^^^ = 
7x3600x(750kms"') - 1.89x10^ km, where we have assumed 
that magnetic cloud moves with the maximum in-situ speed of 
the ambient solar wind behind the shock Vmax - 750 km s"'. 
The magnetic field at the centre of the cloud is estimated to be 
(Wang et al. 2006) Buc = 50 nT . The most likely solar origin 
of this magnetic cloud is a full halo CME that was first recorded 
at 08:50 UT on 2003 November 18 at a height of Rq = 6.3Re = 
2.76x10'^ cm (http://cdaw.gsfc.nasa.gov/CMEJist/). The plane- 
of-sky expansion speed Vexp of the halo CME as recorded in 
the CME catalogue is Vexp = 1660 km s"'. Following Schwenn 
et al. (2005), we computed the radial speed V - 0.88 Vexp = 
1460km s^'. The total time elapsed between 08:50 UT, 2003 
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Table 3. Derived parameters for Forbush decrease for 2003 November 20 event 



Quantity 


NW 


N 


NE 


W 


V 


E 


SW 


S 


SE 


FD magnitude 


0.90 


0.91 


0.89 


1.26 


1.24 


1.04 


1.24 


1.24 


0.88 


FD onset 


13:12* 


11:46* 


02:10* 


01:26 


05:02 


21:07 


06:29 


11:02 


21:07 


FD end' 


04:57 


04:26 


02:32 


06:00 


04:54 


05:00 


04:16 


04:12 


02:54 


cr^ (Kolmogorov) 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 



* These times refer to 2003 November 21. 

' The times for FD end refer to 2003 November 24. 
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Fig. 11. The muon flux in each of the nine directions is shown 
for the Forbush decrease on 2003 November 20. The Unestyles 
and markings are the same as those for fig.|9] 



November 18 (the time the CME was observed at 6.3Rq) and 
10.1 UT, 2003 November 20 (the time the mag cloud was in- 
tercepted by ACE, Wang et al. 2006) is 1.76 x 10^ s. Since the 
ACE satellite is at a distance of ~ 1.5 x 10^ km from the Earth, 
we estimate that it takes ~ 2000 s for the magnetic cloud to ar- 
rive at the Earth at the solar wind velocity V^ax - 750 km s"^ 
Therefore, T = 1.76 x lO'^ + 2000 = 1.78 x 10^ s . 

6.2. Forbush decrease data 

Table 3 gives the data derived from the Forbush decrease ob- 
served on 2003 November 20. This is a comparatively smaller 
event, and there is considerable variation in start times for the 
different bins. This variation is most probably due to uncertainty 
caused by oscillations in the fluxes before the decrease (mainly 
in the NE and N bins), some of which even remain in the fil- 
tered data. In this case since the magnitude of Forbush decrease 
is comparatively small, the diurnal variations are comparable to 
the decrease, and it is not possible to fit the unfiltered data to ob- 
tain any meaningful estimates for the parameters of the Forbush 
decrease; hence, it is necessary to only use the filtered data for 
this purpose. Furthermore, it is very likely that there is some 
CME -related anisotropy remaining in the data, and removing it 
requires the technique used by Munakata et al. (2003, 2005) and 
Kuwabara et al. (2004), which is outside the scope of the current 
paper 



The shock driven by the CME arrived at the ACE space- 
craft at 07:27 UT (Huttunen et al. 2005). Since the average 
shock speed from the sun to the Earth is 1020 km s"' (Wang et 
al. 2006), it would take around 1470 s to traverse the ~ 1.5 mil- 
lion km between the ACE satellite and the Earth, and it would 
have reached the Earth at around 07:51 UT, 2003 November 20. 
On the other hand, the ACE spacecraft first intercepts the mag- 
netic cloud at 10:06 UT (Wang et al. 2006), and the cloud would 
have reached the Earth at around 10:36 UT on 2003 November 
20, taking some 2000 s to traverse the distance between ACE 
and the Earth at the solar wind velocity of V^ax = 750 km s"'. 
Comparing these times with the start times of the Forbush de- 
crease given in Table 3, it is evident that the magnetic cloud start 
time is before the start time in most of the bins, and after it in 
a few of them. However, the end of the Forbush decrease takes 
place well after the entry of the magnetic cloud. Thus the overall 
decrease should have contributions from both the shock and the 
ejecta/magnetic cloud, and we assume that the CME/magnetic 
cloud-related component of the decrease accounts for 50% of 
the total decrease, as discussed earlier in § 3. The estimate of 
the turbulence level cr^ (Eq..[T2ll given in Table 3 was computed 
according to the procedure discussed in § 3. 



7. Results 

We use the observational parameters for each event summarised 
in Tables 1-3 using the procedure outlined in § 3.2 to constrain 
the turbulence level cr^ (Eq. [T2] | for each event. This quantity 
represents the ratio of the energy density in the turbulent mag- 
netic fields to that in the mean magnetic field. We computed cr^ 
separately for each directional bins, and the results are Usted in 
the last row of Tables 1-3. The average value of cr^ for each 
event is quoted in Table 4. 

We note that the turbulence levels near the CME/ejecta front 
for the 2001 April 11 and 2003 November 20 events are fairly 
small (around 2% and 6% respectively). These numbers are 
commensurate with what would be expected in the general qui- 
escent solar wind at a few tens of Rq (e.g., Spangler 2002). The 
turbulence level near the CME/ejecta front of 2003 October 29, 
on the other hand, is significantly higher (around 249%) and is 
representative of fairly strong turbulence. This is not surprising, 
since this event was a much stronger one, and was associated 
with a significantly faster shock. It is conceivable that strong tur- 
bulence associated with the shock front can affect the vicinity of 
the CME/ejecta behind it. 

Our final result regarding the turbulence level cr^ is not very 
sensitive to the precise value of a (Eq.|9]l. A change of a factor 
of 2 in a changes cr^ by 2 to 4%. We also pointed out earlier that 
the estimate for the initial height of the full halo CMEs is made 
from the first observation of the lateral extent of the CME, Rq. 
This is the only concrete observable quantity for halo CMEs, and 
we have verified that an error of 1 Rq in Rq will result in an error 
of around 2% in the estimate of the turbulence level cr^. The de- 
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Table 4. Summary of various observationally derived parameters for the three events we consider, leading to cP- (Eq.[T2li as our final 
result. 



Event 


^0 


^exp 


T 


R(T) 


Bmc 


Average cr'^ (Kolmogorov) 






(kms-') 


(hr) 


(AU) 


(nT) 




2001 April 11 


2.84 


2411 


47.1 


0.10 


17.8 


2% 


2003 October 29 


5.84 


2459 


24.8 


0.14 


44 


249% 


2003 November 20 


6.3 


1660 


49.4 


0.12 


50 


6% 



crease magnitudes vary somewhat from one module to the next, 
as noted earlier. We have verified that these differences in the 
estimate of the Forbush decrease magnitude contribute a scatter 
of only around 2% to our final result for the turbulence level cr^. 
The parameters that influence our final result most strongly are 
the average propagation speed V and the total time T elapsed 
between the time when the CME was first observed and when it 
was intercepted as a magnetic cloud by the ACE spacecraft. 



8. Summary 

We have focused on the contribution of the CME/ejecta in pro- 
ducing Forbush decreases observable from the Earth, using data 
for three events observed with the GRAPES-3 tracking muon 
telescope at Ooty. We selected events that clearly have their ori- 
gins in full-halo CMEs that originate close to the centre of the 
solar disc. The Forbush decreases produced by these events will 
therefore be due to the shock, as well as to the CME/ejecta. 
While it is somewhat difficult to precisely separate out the shock- 
associated and CME/ejecta-associated contributions, it is rea- 
sonable to assume that the CME/ejecta effects accounts for 50% 
of the overall Forbush decrease. We used CME measurements 
from LASCO and magnetic cloud measurements from WIND 
and ACE to obtain various geometrical parameters of the CME 
and its near-Earth manifestation. The CME is a substantially 
closed magnetic structure that starts out with a negligible num- 
ber of high-energy protons inside it. We consider the diffusion 
of high-energy cosmic -ray protons into the CME as it propa- 
gates and expands on its way from the Sun to the Earth. Since 
the diffusion primarily takes place across the magnetic fields that 
enclose the CME, we use an expression from the literature for 
the cross-field diffusion coefficient that is derived from exten- 
sive simulations of cosmic rays propagating through turbulent 
magnetic fields. The diffusion coefficient depends upon several 
quantities like the ratio of the proton Larmor radius to the max- 
imum scale length of the turbulence p (Eq. ifTTI ) and the level 
of turbulence cr^ (Eq. |[T2l ). The turbulence near the CME is ex- 
pected to result in a magnetic-field line random walk, leading to 
scattering of protons and contributing to their cross-field diffu- 
sion. 

We constrained the magnitude of cr^ (which is the ratio of the 
energy density in the random magnetic fields to that in the large- 
scale magnetic field) in the vicinity of the CME front for each of 
the events we analysed. We find that the turbulence level o"^ near 
the CME front for the 2001 April 1 1 Forbush decrease event is 
~ 2% and ~ 6% for the event of 2003 November 20, while it is 
~ 249% for the much more energetic event of 2003 October 29. 
These estimates may be regarded as an average over the jour- 
ney of the CME from the Sun to the Earth, for we do not take its 
possible radial evolution into account. The radial evolution of cr^ 
is unlikely to be appreciable, at least in the outer corona, as we 
discuss below. If we take turbulent magnetic field fluctuations 
to be representative of those in the electron density (Spangler 
2002), we note that the power-law index for electron density 



fluctuations AA^g as a function of heliocentric distance ranges 
from -2.2 (Manoharan 1993) to -1.7 (Faflows et al. 2002). This 
is fairly close to the power law index of the background electron 
density A^^,, which is around -2 in the outer corona (e.g., Leblanc 
et al. 1998). In other words, the quantity ANe/Ne does not vary 
appreciably with distance from the Sun, and we can surmise that 
cr^ behaves similarly. 

To the best of our knowledge, these are the first quantitative 
estimates of the turbulence levels near CME fronts. Manoharan 
et al. (2000) have noted that the turbulent density spectra of 
the plasma near CMEs are significantly flatter than what is ob- 
served in the slow and fast solar wind, and are similar to what 
exists in the near-Sun solar wind acceleration region. However, 
they were unable to comment on the level of the turbulence. 
Badruddin (2002) suggests that the level of turbulence near 
the CME/magnetic cloud might influence the magnitude of the 
Forbush decrease. However, he only considered the local effect 
of the turbulence, while our treatment takes into account the cu- 
mulative effect of the turbulence near the CME on the cross- 
field diffusion of protons into it. Wibberenz et al. (1998) quote 
the results of Vanhoefer (1996) who makes a simple estimate of 
the cross-field diffusion coefficient by calculating the "filling in" 
timescale of the CME by energetic protons. This method results 
in a rather low estimate for the diffusion coefficient, and they 
speculate that this might be due to the smooth magnetic fields 
enclosing the CME. Our treatment considers this aspect in de- 
tail, incorporating the expansion of the CME, as well as results 
for the cross-field diffusion coefficient adopted from extensive 
numerical simulations. 
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